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Lepton flavor violation (LFV) has been observed in neutrino oscillations. For charged lepton FV 
decays only upper limits are known, but sizable branching ratios are expected in many neutrino 
mass models. High scale models, such as the classical supersymmetric seesaw, usually predict that 
decays k — >■ 3/^ are roughly a factor a mailer than the corresponding decays k — > lj"f. Here we 
demonstrate that the Z^-penguin diagram can give an enhancement for decays k — >■ 3lj in many 
extensions of the MSSM. We first discuss why the Z°-penguin is not dominant in the MSSM with 
seesaw and show that much larger contributions from the Z^-penguin are expected in general. We 
then demonstrate the effect numerically in two example models, namely, the supersymmetric inverse 
seesaw and R-parity violating supersymmetry. 
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Introduction: Neutrino oscillation experiments [TJ |2] 
have firmly established that lepton flavor is violated in 
the neutrino sector, with two of the three measurable 
mixing angles being surprisingly large. Observation of 
the characteristic "neutrino dip" leaves no doubt that 
neutrinos have mass [31 H] and quite accurate values for 
the mass squared differences are known now [5 . In the 
charged lepton sector, however, only upper limits on LFV 
branching ratios, such as /i — 67 [16 or /i 3e [7 , exist. 

Extending the standard model (SM) only by neu- 
trino masses does not automatically lead to measurable 
charged LFV (CLFV), but sizable branching ratios are 
expected in many models. In fact, on quite general 
grounds one expects large CLFV, if physics beyond the 
SM exists at the TeV scale. A prime example for this 
observation is supersymmetry (SUSY). Here, the mass 
matrices of the new scalar particles need not (and in gen- 
eral will not) be aligned with those of the SM fermions. 
CLFV decays will occur and one can estimate roughly 
the branching ratios for radiative decays as |8j 
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where {m?^)ij parameterizes the dominant off-diagonal el- 
ements of the soft SUSY breaking slepton mass matrices 
and rnsuSY is the typical mass of the SUSY particles, 
expected to be in the ballpark of 0(0.1 — 1) TeV. Rather 
small off-diagonal elements are required to satisfy exper- 
imental bounds |7] . 

In the Constrained Minimal Supersymmetric exten- 
sion of the SM (CMSSM), on the other hand, CLFV is 
zero, just as in the SM, simply because neutrinos are as- 
sumed to be massless. Extending the CMSSM to include 
neutrino masses (and mixings) , for example by a seesaw 
mechanism, then leads to CLFV decays, because the fla- 
vor violation necessarily present in the Yukawa couplings 
is transmitted to the slepton mass matrices in the RGE 



(renormalization group equation) running \W . In such 
high-scale neutrino mass models, with only MSSM par- 
ticle content at the electroweak scale, it has been shown 
that the photonic penguin diagram gives the dominant 
contribution to li — )• 3lj decays in large regions of pa- 



rameter space, 
derived [TT] 



In this case a simple relation can be 
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Thus, usually it is concluded that the decays li Ij'^ are 
more constraining than the decays U — >■ 31 j. 

Apart from the photonic penguin, there are also box 
diagrams, Higgs- and Z°-penguin contributing to the de- 
cays li — >■ 31 j. The latter diagram is not per se smaller 
than the photonic penguin. Rather, as we will show be- 
low, in models with only the MSSM particle content and 
couplings, the Z'^-penguin is suppressed by a subtle can- 
cellation among different terms in the amplitude. Such 
a cancellation, however, can be easily spoiled if there 
are (a) new couplings and/or (b) a larger particle con- 
tent than in the MSSM. Then, as we will discuss, the 
.Z'^-penguin can easily give the dominant contribution to 
k — )■ 31 j. We will demonstrate this fact numerically with 
two typical example models: (i) a supersymmetric inverse 
seesaw and (ii) R-parity violating SUSY. The former is 
an example of a model with extended particle content, 
while the latter is an example of a model with the MSSM 
particle content but new interactions. As we will show, 
in such models li — >■ 3lj can be more constraining than 
h IjJ- This is the main result of the present paper. 

Analytical discussion: The total width of the li — >■ 3lj 
decay contains contributions from the photon penguin. 



^ An exception from this rule is the decay t 
large tan ^ [TO) . 



3fj, in the limit of 
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the Higgs penguin, the Z°-penguin and boxes. Consid- 
ering only contributions from photon and Z'^-penguin. 
which are the ones of interest to us and numerically 



the most important ones, the total width T = r{l- ^• 
Ijijlj^) can be written as |12| : 
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{2 {\Fll? + \Frb\^) + \Flb\^ + \Frl\^]+Iaf 
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Here, terms denoted A (F) are due to photon (Z°) ex- 
change and Iaf denotes their interference terms, irrele- 
vant for the following discussion. Both, photon and Z° 
penguins have chargino-sneutrino and neutralino-slepton 
contributions. Exact definitions can be found in [11 . We 
will focus on the chargino loops for brevity here, since the 
effects we are interested in are most pronounced in these 
loops. The photon contributions are 

^(e)L.H ^ -Lo^f .(.^) (4) 

m- 

whereas the Z-contributions read 

Fx = -^^^ y'^f'^V) (5) 

with X = {LL, LR, RL, RR}. In these expressions Oy'^ 
denote combinations of rotation matrices and coupling 
constants and s(x^) and t{x'^) are short-hands for the 
Passarino-Veltman loop functions which depend on = 
rn^^/rn^. For precise definitions see 

The scaling A ~ mg^gy and F ~ can be un- 

derstood, in principle, from simple dimensional analysis. 
The width of the decay is proportional to mf, , so both A 
and F must he A, F oc m^^. In penguin diagrams it is 
the smallest mass term in the loop which sets the scale, 
which in F is mz ■ Due to the masslessness of the photon 
in case of A the smallest mass scale in the loop is msusY ■ 
With TO^ ''^'susY penguin can, in principle, be 

even more important than the photonic one. 

Numerically, however, it has been found in case of the 
MSSM that the photonic penguin is dominant pTj. This 
can be understood as follows. To simplify the discussion, 
we neglect first Fr, since it is always proportional to the 
charged lepton Yukawa couplings. Consider then n gen- 
erations of sneutrinos and neglect the chargino mixing. 
In this simplified scenario only the wino contributes to 
F^'-'' and it can be written as 

i^i'^ = A/wavc + Mpl + Mp2 (6) 



with 

^wave ^5i(52CW -5lSw)2'y ^yVwave (7) 

Afpi = -glcwZ\^Z^*f;, (8) 
Mp2 - l9l{g2Cw+9iSw)Zl^Z'^*f;^ (9) 

Summing over the index i is implied. The terms 
in the sum come from different types of diagrams: 
wave function diagrams (Mwavc), penguins with the Z'^- 
boson attached to the chargino line (Mpi) or the sneu- 
trino line (Afp2). Moreover, Cw = cosOw, sw = 
sin^VF- Zy is a n x n unitary matrix that diag- 
onalizes the mass matrix of the sneutrinos and we 
used the abbreviations /^avc = ~^i('^|± i "^1;); /pi = 
5C'o(TO?^,m|±,m?±) - m^±Co(TO?^,m|±,m|±),/;2 = 
iC'o(TO|±,rn,?.,m?.). The sum in eq. ^ vanishes exactly 
as can be seen by grouping the different terms 

F^'^ = lgicwZ\,'Z'^*Xl + \glg^swZ^^ZfXl (10) 

with X\ = - 2/^1 + f;^,X\ = fl^ - Using 

the exact expressions for the loop functions [TT] one finds 
that the masses cancel out and these combinations be- 
come just numerical constants: XI — —j and X2 = —\. 
Therefore, one is left with fI^""^ oc J2i ZyZ'^*, which van- 
ishes due to unitarity of the Zy matrix. 

This cancellation can be spoiled by two effects, either 
(1) the sneutrinos mix with other particles which are 
not SU{2)l doublets so that the factorization no longer 
holds, or (2) the charginos are not pure wino and higgsino 
states. The last effect is of course present in the MSSM 
and therefore this cancellation is not exact. Nevertheless, 
the ^''-contributions are suppressed due to their propor- 
tionality to the square of the chargino mixing angle (two 
wino-higgsino insertions are necessary since there is no 
— vl — II coupling). We neglected so far Higgsino 
interactions because in many models these couplings are 
very small in comparison to the gauge interactions (for 
example, a SUSY scale type-I seesaw model would have 
^ 10~^). However, in models where the Higgsino can 
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have much larger Yukawa interactions, a large enhance- 
ment of the Z'^-contributions can be expected. This will 
be addressed numerically in the next section. 

Before turning to a numerical discussion, consider for 
simplicity a toy model consisting of two generations of 
left-handed sneutrinos which can mix with one genera- 
tion of right-handed sneutrinos. A 3 x 3 rotation matrix 
is in general parametrized by 3 angles, but we will as- 
sume here for simplification that two of them vanish and 
call the third one In addition, we introduce a new 
interaction for the Higgsinos kv'^HuIl- We give in Fig. 1 

(c) 

the computed F£ for arbitrarily chosen sneutrino and 
chargino masses as a function of ^ for different values 
of K. The dashed line shows the case for = k = 0. 

(c) 

As clearly seen, F£ depends on the left-right mixing al- 
ready for small values of k. However, increasing k, F^^ 
becomes totally dominated by the new k interactions and 
enhances Br(/i — )■ 3lj). 
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FIG. 1: F^''' for our toy model as a function of the sneutrino 
left-right mixing angle ^ and for different values of k: 10"'' 
(blue), 10"^ (red), 0.1 (orange) and 1.0 (black). 

Numerical examples: We turn to the full fledged nu- 
merical study of two examples: an inverse seesaw model 
and the MSSM with i?-parity violation. For this pur- 
pose, we have created for both models SPheno modules 
[131111] using the Mathematica package SARAH [rSHTf]. 
These modules calculate the low-energy observables ex- 
actly including all possible diagrams [17J. 

Inverse Seesaw: In inverse seesaw the MSSM particle 
content is extended by three generations of right handed 
sneutrino superfields and of gauge singlets Ns which 
carry lepton number |18| [T9] . The superpotential reads 

Wis = WMSSM+Y^i^^LH^+MnP-Ns + ^NsNs (11) 

After electroweak symmetry breaking (EWSB) the effec- 
tive mass matrix for the light neutrinos is approximately 
TOy ~ ^Y^{M]^)~^ ^jqM~^Yj . Since /za? can be of 
0(10"^) keV or even smaller while Mn is of 0{msusY), 
the neutrino Yukawa couplings have to be much larger 
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FIG. 2: Top: Different contributions to Br(/i — > 3e) as 
function of the lightest sneutrino mass: Z^-penguins (red), 
photonic penguins (blue), combined Higgs-penguins/box dia- 
grams (green). Bottom: Br(/x — >■ 3e) (black) and Br(/i — > 67) 
(blue) and the current experimental bounds (dashed lines). 
The dips are an effect of a mass crossing between charginos 
and sneutrinos. 

than for a standard weak-scale seesaw to explain neu- 
trino data. 

Due to the extended particle content, new contribu- 
tions for Br(Z.i — > 3/j) are expected in the inverse seesaw. 
For example, the Higgs mediated contributions were re- 
cently studied in [20:. In Fig. |2] (top) we show the dif- 
ferent contributions to Br(/x — >■ 3e) for a variation of the 
SUSY masses. To disentangle RGE effects we have calcu- 
lated once the spectrum for a GUT input (mo — 500 GeV, 
Ml/2 = 1 TeV, tan(^) = 10, Aq = -300 GeV) and 
rescaled all dimensionful parameters at the SUSY scale. 
This changes the sfermion masses but not the mixing 
matrices. has been chosen to explain neutrino data 
for diag(/iAr) = 10"! keV and Mr = I TeV. Clearly, 
the Z°-penguins dominate and are nearly independent of 
the SUSY scale. Only in the limit msusY ~^ fnz the 
other contributions can compete. In Fig. [2] (bottom) the 
branching ratios for /i — > 67 and /i — > 3e and the current 
experimental bounds of 2.4 • 10~^^ and 1.0 • 10~^^ are de- 
picted ^iZj. While Br(/i — > 67) would be in conflict with 
experiment only for m^j^ < 1.2 TeV, Br(/i — > 3e) rules out 
the entire range. In this example, we have assumed 
and Mr to be diagonal and all flavor violation comes from 
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FIG. 3: Br(^ — > 3e) (black) and Br(/i — > e-y) (blue) varying 
A132 ■ A232 in a mixed bi- and trilinear (solid lines) and a pure 
tri-linear (dotted lines) RpV scenario. The dashed lines show 
the experimental limits. 



Yy, as is usually done in literature. However, neutrino os- 
cillation data could equally well be fitted with the flavor 
violation coming from fi^ and Mj^. In that case CLFV 
observables would be much smaller and consistent with 
experimental data. However, the relative order between 
2- and 3-body decays won't change, i.e. Br{li — )• 3lj) 
will be most likely observed before Br(^i — > /j7) if inverse 
seesaw is realized in nature. 

R-parity violation: As second example, we take the 
MSSM particle content but extend the superpotential by 
the lepton number violating terms \22\ 



1 



1 



= Wmssm + -KjkULkEl + -A, 



(12) 

While the e-parameters are highly constrained by neu- 
trino data |21| . the bounds for the tri-linear couplings 
are much weaker and some entries can be of 0(1) |23| . 
In the following, all entries of A and A are set to zero 
but A132 and A232. We give in Fig. [3] the results for 
Br(/i — )■ 3e) and Br(/x ej) for a mixed bi- and tri- 
linear as well as for the pure tri-linear scenario varying 
A132 ■ A232- In the mixed case and v]^ have been cho- 
sen to be consistent with neutrino data and a moderate 
flavor violation in the sneutrino sector has been induced 
by m? ^ = (45 GeV)^. It can be seen that in the mixed 
case Br(/i — )■ 3e) > Br(/i — > ej) holds when A132 ■ A232 
crosses 2.5-10^^, while for the pure tri-linear case without 
any flavor violation at tree level in the sneutrino sector 
the three body decays dominate even for much smaller 
values. 

In both cases we get an upper limit for A132 ■ A232 of 
2.5 • 10^^ from the bounds on Br(/x 3e) for sneutrino 
masses of 730 GeV. So far, in the literature just the limits 
for = 100 GeV from the photonic penguins |24j have 
been published. These are much weaker, after rescaling 



penguin can give the dominant contribution in lepton 
flavor violating three body decays in many models. The 
importance of the Z°-penguin increases with increasing 
SUSY particles masses. As numerical examples, we have 
briefly discussed the supersymmetric inverse seesaw and 
the MSSM with i?-parity violation. 
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> 100 GeV; 

Summary: We have shown in this letter that the Z'^- 
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